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ABSTRACT 

We have probed the tidal tails of the Sagittarius dwarf spheroidal galaxy using 
an as yet unmeasured sample, stars in the Red Clump. While contamination from 
the disk of our Galaxy is a primary concern and a possible reason for ignoring this 
population of stars until this point, we have developed a method to systematically 
eliminate this background using the stellar parameters of our program stars. We 
present this procedure for selecting Red Clump stars. Using this method, we 
measure the mean of the metallicity distribution of the streams in two areas to 
be [Fe/H] = -0.4 and -0.2 which agree well with previous measurements of the 
core of the dwarf galaxy. We also find the kinematics of one sample to have 
a velocity of -120 ± 20 km s^^ which is in good agreement with an M giant 
measurement in the same area. Finally, we find that the density of Red Clump 
stars varies as a function of distance from the Sagittarius dwarf galaxy. Densities 
of p = 3.0 ± 0.3 RC stars/kpc^ and p < 1.3 ± 0.3 RC stars/kpc^ in areas at 
different distances from the core of the dwarf galaxy were found. This factor 
of two difference is likely the result of the Red Clump stars being preferentially 
stripped during the most recent peri-galacticon passage. This suggests that the 
Red Clump stars may have been more centrally located in the Sagittarius dwarf 
galaxy than their Population II counterparts. Observations were carried out at 
Kitt Peak National Observatory using the WIY J3 telescope. 

Subject headings: galaxy: halo — galaxies: interactions — stars: horizontal branch 



^The WIYN Observatory is a joint facility of the University of Wisconsin-Madison, Indi- 
ana University, Yale University, and the National Optical Astronomy Observatories. 
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Introduction 



Since its discovery ( llbata et al.lll994| ). the Sagittarius dwarf spheroidal galaxy (Sgr) 
has served as an important probe for the understanding of the properties of our Galaxy. 
Unlike other previously identified Milky Way satellite galaxies, Sgr is a relatively nearby 
galaxy with clear evidence of being tid ally disrupted. Less than two years after its discovery 



the Sgr was found to have 'tidal tails' (IMateo et al. 



19961 ) - stars which have been stripped 



from the main body of the dwarf galaxy by its interaction with the Milky Way Galaxy. 
This led to the finding that the Sgr actually ha d an extended and pro minent stream of stars 



that was ultimately traced over the entire sky (IMajewski et al 



20031 ) 



The core of the Sgr is located below the bulge of our Galaxy as seen from the location 
of the Sun. Its orbit has been found to be near polar and both leading and trailing streams 
of tidally stripped stars accompany the core. These streams trace out the full orbit of the 
Sgr and overlap each other in certain places. Unknown at this time is the number of orbits 
the streams represent and how long ago these streams have been pulled from their host. The 
Sgr and its extended tidal streams of stars are a very important and unique opportunity to 
probe several areas of modern astrophysics. By studying the locations and kinematics of 
the tidal streams, one is probing the orbit of the dwarf galaxy as it moves in and around 
our Galaxy's gravitational potential. The orbit of such a system is quite sensitive to the 
physical properties of the potential and more importantly on the distribution of the mass 
that creates the potential. Precise mapping of the orbit of the Sgr will therefore lead to a 
constraint on the shape of the dark matter that surrounds our Galaxy. 



A photometric study by 



Bellazzini et al. 



( 120061 ) found a gradient in the ratio of old 



to young stars along the arms of the Sgr tidal streams. In their work, they looked at the 
number of Blue Horizontal Branch (BHB) stars compared to the number of Red Clump 
(RC) stars in the core of the Sgr and in an area far from the core along one of the tidal 
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streams. They found that the BHB stars are five times more abundant in the stream 
than in the core of the dwarf galaxy when comparing to the RC stars. Their results 
were based solely on star counts from photometric data which means contamination from 
various sources (galaxies in the BHB colors and foreground disk stars in the RC colors) 
were harder to eliminate than if there were spectroscopic information available. Seeing this 
as an opportunity, we set out to observe RC stars spectroscopically and try to confirm 
their membership in the tidal streams of the Sgr kinematically. The benefits of these 
observations would be two- fold. First, it would give a more precise measure of the number 
of intermediate age stars in the streams of the Sgr. Second, the kinematic measurements of 
the stream would be directly comparable to previous measurements using a new and unique 
stellar population and could help constrain different models of our Galaxy's dark matter 
halo. 

In section [2] of this paper we will discuss the sample selection criteria and the 
observations and reductions of our spectroscopic data. Section [3] describes the analysis 
procedure used to provide a clean RC sample. In section H] we present the results found 
using this stellar sample and follow it up in section |5] with a discussion of these results. 



Data 



2.1. The Sample 



The candidate stars were selected usin g the Sloan Digital Sky Sur vey (SDSS) online 



database 



Using 



previo us measurements ( iMajewski et al. 



2003 



models ( jLaw et al. 



20041 ) and simulated 



20051 ) of the Sgr tidal tails, two locations for study were chosen. These 



points are dynamically very similar according to the models but are at different distances 



2 



http : //www . sdss . org 
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from the Sgr. The models predict the possibihty of more than one arm of the tidal debris 
overlapping in these areas that are at roughly the same distance from the Sun but have 
velocities which are different by several hundred kilometers per second. This point is 
important for several reasons. The most important reason is that only a single magnitude 
range would be necessary to isolate stream stars from different arms. Also, the velocity 
difference is necessary to separate one tidal arm from the other. Furthermore, a large 
difference in velocities, with one being positive and the other negative, would help to ensure 
a clean sample. Table [1] shows relevant values for our selected fields and Figure [1] has 
our field position s in Galactic latitude and longitude overlaid on the spherical model of 



Law et al. 



(120051 ). 



To choose candidate RC stars, color and magnitude cuts were determined that would 
minimize the amo unt of contamination f r om va rious sources (mainly the disk of our Galaxy) . 



Using the work of 



Lay den fc Sarajedinil (|2000| ) we were able to obtain a reasonable color 



cut for our sample. Their photometric data of the globular cluster M54 and Sgr field stars 
provided a direct measure of the RC color for the Sgr proper which is given as V — I = 1.13 
in their paper. To transfer to the SDSS photometric system we used the photometric 
equations by Robert Lupton given on the SDSS websit f|. The magnitude range was 
influenced by the expected distances of stream stars, telescope aperture, and the reliable 
range of the SDSS. The selection criteria adopted for this research was 0.5 < {g — r)o < 0.65 
and 15.5 < Qo < 18.7. A go versus {g — r)o CMD of stars in our fields with our selection box 
overlaid can be found in Figure O The range in magnitude allows for two things. One is 
the possibility there is more than one wrap of the tidal debris and the wraps are at slightly 
different distances. The other is the possibility suggested in the models that there is a 
portion of the stream at a much closer distance to the Sun. 



http : //www . sdss . org/ dr6/ algorithms/ sdssUEVRITraJisf orm . html 
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The chosen color and magnitude ranges will undoubtably make contamination from the 
disk of the Galaxy an important concern. The elimination of this background is possible, 
however, if the stellar parameters and kinematics are known. The disk stars will have a 
different velocity signature than the Sgr tidal stream stars will. The disk stars should also 
have higher surface gravities than the stream stars since they should still be on the main 
sequence. Combining this information together will give excellent discrimination between 
the background disk stars and the desired stream stars as will be shown. 

2.2. The Observations 

Proposals for the observations of the candidate stars chosen as described in section 
12.11 were submitted to the NOAO. This research was granted time two semesters in a row 
using the WIYN telescope at KPNO with the Hydra instrument. Observations were carried 
out in August 2007 and February 2008. The telescope and instrument combination was 
critical since the necessary statistics can be obtained in a much shorter timespan because of 
the ability to take many stellar spectra in the same integration. For our research the blue 
fibers and Simmons Camera were used with no filter and a 600 lines per mm grating. This 
gave a resolution of approximately 4.6 A and a dispersion of just over 1 A per pixel. Our 
detector setup provided us with a wavelength coverage of roughly 4000-5500 A. The density 
of candidates and the positioning constraints of the instrument allowed between 43 and 62 
objects to be observed in a given field. 

On each program field we took four integrations of 30 minutes each and co-added them. 
This gave us the two hour integration times needed to get a decent signal on our faintest 
stars and helped eliminate difficulties with cosmic rays. We found that with this program 
we got typical S/N ratios of ~ 45 at a go magnitude of 16.5 and ~ 20 at a go magnitude of 
18. 
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The reduction of all images was performed using the IRAF software packagcl. Bias 
levels, CCD response, fiber locations, and dispersion were all taken into account in the 
dohydra routine of IRAF. This facilitated the extraction of individual spectra from the 
images taken at the telescope. 



3. Analysis 

In an effort to extract meaningful information from the spectral data obtained as 
described in section [21 an analysis software framework was developed by the authors. 
The ability to separate the RC stream stars from main sequence disk stars hinged on the 
ability to determine the kinematics of our sample as well as reliable stellar parameters. Of 
particular interest is the ability to quantify surface gravity since this parameter is a key 
discriminate for the RC sample. 



3.1. The Procedure 



Our procedure uses a grid of synthetic spectra to compare to our observed data. The 
synthetic stellar spectra were generated using the SPECTRUM prograno with ATLAS9 
model atmosphereqj as input. These grids were generated to emcompass the entire 
parameter space likely to be observed in the program stars. The grid size for the effective 
temperature was 4000 to 6000 K with a step size of 250 K. For surface gravity a step size of 
0.5 dex was used over a range of 0.5 to 5.0 dex. The metallicity range was -2.5 to +0.5 dex 



"^jhttp : / /iraf . noao . edu| 

^Developed by Richard Gray, http : / / www . phys . appstate . edu/spectrum/ spectrum . html 
^From Robert Kurucz, http : / /kurucz . harvard . edu 
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in increments of 0.5 dex. 

A primary concern we had with the synthetic data (as well as the scientific data) 
was how the normalization would be carried out. SPECTRUM allows for the synthetic 
spectra to be output in a normalized form. This true normalization can be quite different 
from the normalized, pseudo-continuum found for the science data by fitting the non-fiux 
calibrated spectra. Direct comparison between the two normalizatons leaves open the 
possibility that systematic errors could be introduced in the form of incorrect line strengths 
because of incorrect continuum placement. To avoid this possible problem we decided to 
output the synthetic spectra in absolute fiux values and normalize them later in the same 
way we normalized the scientific data. The normalization procedure we used involved the 
ROOT software packagqj. There is a routine called Background in the TSpectrum class 
of this package that can determine a background level of a spectrum. To determine the 
background, the routine uses a Sensitive Nonlinear Iterative Peak clipping algorithm. More 
information about the routine and the algorithm can be found in the reference material on 
the ROOT webpage. Normalizing both the synthetic and scientific spectra with this routine 
meant that any incorrect placement of the continuum would be in both sets of data and 
therefore would not lead to systematic errors when comparing the two. 

The method determined to be best suited for determining the stellar parameters of our 
program stars is a type technique. The difference between a scientific and synthetic 
spectrum was calculated for a specific set of synthetic stellar parameters. That difference 
was then compared to the difference from other sets of parameters. A convergence criteria 
was adopted to rehably determine the best parameters that employed the method of 
steepest-descent. Beginning with expected median values for our sample, three parameters 
were allowed to vary sequentially and independently to locate a new minimum in the 



ROOT is a data processing framework available for download at http : / /root . cern . ch 
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difference of the neighboring points of the grid space. Convergence was accomphshed and 
the parameters considered correct when the resulting values matched the input values (a 
minimum was found). If the resulting values did not match the input values the process 
was repeated using the new resulting values as input values until convergence was reached. 
Once the best grid point was found, the better result from a cubic spline and quadratic 
fit was used to determine the most accurate solution from the minimum value and its 
neighboring values. 

Since the full spectrum of a star is not sensitive to all stellar parameters at every 
single wavelength, specific sections of the full spectrum were selected to help determine 
each parameter. These different areas were given weights in order to include as much useul 
information as possible without introducing too much noise. A table of the wavelength 
ranges, adopted weights, and constrained stellar parameter can be found in Table [H One 
sees that we are using the Balmer lines as an effective temperature indicator with greatest 
weight on the line. We have done this because the S/N drops rapidly at the blue 
end of the spectrum limiting the usefulness of the and H5 lines. Surface gravity was 
determined from the wavelength region that included the magnesium B-lines and the 
magnesium- hydride molecule. Other areas were investigated but nothing was found to 
improve the results from this single range. Our abundance information is spread over the 
entire wavelength range but excludes the Balmer lines and magnesium B-line regions. All 
the areas between our other indicators were used with greater weight on the red end of the 
spectrum because of the higher S/N and prominent iron lines in this region. 



3.2. The Validation 



To confirm that our software produc e d qua lity results, a sample of approximately 70 
stars was used from lAUende Prieto et al.l (|2008[ ). These stars have spectra obtained from 
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the SDSS telescope and stellar parameters determined with the pipeline software developed 
for SEGUE. They have also been observed at high resolution and stellar parameters have 
been extracted from these results as well. The parameter range of these standard stars 
covers the RC region we are interested in. The effective temperature of the sample ranges 
from roughly 4500 to 6250 K. The surface gravity ranges from roughly 1.0 to 5.0 dex 
and the metallicity ranges from roughly -2.5 to 0.0 dex. As can be seen in Figure [31 the 
procedure developed for this research matches the high resolution results quite well (not 
shown is the comparison between our results and the SEGUE pipeline results which shows a 
similar agreement). The figure shows the difference between stellar parameters determined 
using our routines on the SDSS spectra and those determined from the high resolution 
spectra. Of note is the fact that the stellar parameters determined with our technique 
show differences that have means of close to zero and distributions that are Gaussian in 
shape. Furthermore, the widths of the distributions are on the same order as the synthetic 
data's grid spacing which means our resolution is dominated primarily by the grid spacing 
as expected. This means our procedure gives reliable stellar parameters and does not 
introduce any significant systematic errors. Figure [3] also shows the difference in the radial 
velocity determination between the two methods. The width of this distribution is quite 
respectable at 3-4 km s~^ and there is a slight offset of ~ +6 km s~^. Table [3] shows how 
our procedure compares with the high resolution and pipeline results. In summary, our 
software is considerably faster than the SEGUE pipeline and produces results with a similar 
accuracy and precision in this stellar parameter range which makes it an excellent tool for 
the anaylsis of our program stars. 



4. Results 



As was described in section 12.11 our selection criteria allowed for the inclusion of 
quite a significant background, namely the disk of our Galaxy. Since we observed stars 
spectroscopically, however, we were able to take advantage of the fact that the stream and 
disk stars are in different stages of the stellar life cycle. This coupled with the velocity 
information provides discrimination between the stream and disk stars. We will first 
describe the procedure used to select the stream stars from our full sample and will follow 
that with the results these stars provide us. 



4.1. Separating the Sample 

Figure H] shows a plot of GSR velocity versus surface gravity for our program stars. 
This plot shows exactly what is expected: a large number of stars with low velocities 
(~0 km s~^) and high surface gravities (> 4 dex in log g) as well as a population with lower 
surface gravities (< 2.5 dex) and very different velocities (~ -125 km s~^). The former 
group is undoubtedly the disk stars in our sample and the latter is the RC stars in the tidal 
stream of the Sgr. One could put a hard cut on these two variables to produce a reasonable 
sample of stream stars. We chose to use other information to select our sample and confirm 
that we are choosing the right group by reproducing this plot with the sample separated. 

As a start we looked at the dependence of surface gravity on temperature. The results 
can be seen in Figure O The disk-like population is apparent in the lower left corner as 
the large group centered at a temperature of roughly 5500 K and a surface gravity of 
roughly 4.25 dex. The RC stars we are interested in should be above this group in the 
plot (lower surface gravities) and so we adopted the two selection cuts represented by the 
lines in the figure. These cuts are 1) a hard cut at log g < 3.25 (the solid line) and 2) the 
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temperature dependent cut log g < temperature / 500 - 7.55 (the dashed hne). The second 
cut was selected to eliminate stars along theoretical isochrones that would be expected for 
a population of stars similar to the disk of our Galaxy. 

Using these two cuts, we then looked to see if there was any dependence of the GSR 
velocity on the Qo magnitude of the stars. Since the RC is a part of the HB, the absolute 
magnitude of all these stars should be roughly the same. Therefore, by looking at the 
apparent go magnitude we are indirectly examining the distance of the stars. One can see 
in Figure [6] there are two distinct groups of stars. The brighter stars have quite a different 
velocity distribution than the fainter stars do. We interpret this to mean that the brighter 
stars are likely thick disk or halo stars in our sample since they are brighter (closer) and 
have velocities consistent with either of those two populations of our Galaxy. This led us 
to another cut of 3) Qo > 17 (the faintest stars in our data set). If the brighter stars are 
a contamination from the thick disk or halo of our galaxy we would expect them to be 
randomly distributed through our data. If the brighter stars are part of a closer stream 
they would likely be more localized and coherent. Figure [7] shows that indeed the former is 
true and gave us confidence that this selection criteria is sound. 

Armed with the three cuts mentioned, we revisit Figure H] but now make the plot with 
two distinct groups from the sample. The result is Figure [HI Those meeting the criteria 
explained above have been plotted as solid squares and those not meeting the criteria have 
been plotted as crosses. There is obviously a little overlap between the two groups but 
this fact is to be expected. There are RC type stars in different populations of our Galaxy 
and they could have velocities that are close to the velocity of stars in the Sgr stream. 
The significant point is that our systematic selection of low surface gravity stars with 
magnitudes between 17 < go < 18.7 gives us a sample of stars with very similar velocities. 
Furthermore, the velocity at which these stars move is very different from any Galactic 
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component in the same direction. 

The plots discussed in this section are for the observing run in August 2007 only. 
Similar plots for our February 2008 observing run have been made and the same set of 
criteria applied. Now that we are convinced that our data does have a sample of stars that 
are consistent with a RC population in the Sgr tidal streams and we can select them using 
stellar parameters obtained from our analysis software, we can perform measurements of 
the streams. 



4.2. Properties of the Sgr RC 

The number of stars in each group of our data provides some meaningful information 
and is found in Table |H As can be seen, our observations probed a little less than half the 
stars SDSS says there are in this region of the sky with these color and magnitude cuts. In 
our observed sample 14 percent of the August 2007 data and 20 percent of the February 
2008 data are stars we identify as being RC stars. The candidate selection and placement 
of stars on fibers of the instrument is more or less completely randono so we assume our 
sample has no selection bias. Therefore, if we were to completely probe these regions of 
the sky, we would expect roughly 13 and 15 RC Sgr stream stars per square degree to be 
present in the August 2007 and February 2008 data sets respectively. 

Two of our three stellar parameters were used to select our sample so the only unbiased 
parameter left to study is the metallicity. Figures |9] and [TOl show the metallicity distribution 
for the August 2007 and February 2008 datasets respectively. These distributions have been 



^Priority can be put on certain stars but that does not guarantee their placement. It only 
sets an order the software uses to attempt to configure the field and makes little practical 
difference. 
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fit with Gaussian functions and we have found a mean of [Fe/H] = -0.4 for the August 2007 
sample and a mean of [Fe/H] = -0.2 for the February 2008 sample. Both have widths on 
the order of a dex. These measurements are in good agreement wi th what has been found 
previously in the core and tidal tails of the Sgr (see the results of IChou et al.l (120071 ) and 
their discussion in Section 2 of the paper). 

The stream stars in the August 2007 sample are a very distinct group kinematically 
that we can measure quite precisely. Normalized histograms of velocities of the two samples 
we observed (stream and disk) as well as a normalized histogram of stars from the Besancon 
stellar population synthesis model of the Galaxjo are shown in Figure [TTl One sees that our 
disk sample matches very well what is expected from the model. In Figure [T2] we look at our 



samp 



e of RC stars between 1 



(120041 ) and 



Yanny et al. 



Majewski et al. 



5° < / < 175° and compare with the results of 
(j2009l ) in the same region. We have measured the RC component of 
the Sgr stream in this area to have a GSR velocity of -120 ± 20 (stat) km s"^ which agrees 
very well with the M giant measurement of -130 ± 20 (stat) km s^^. The BHB results of 



Yanny et al. 



(120091 ) possibly have a population other than those in the Sgr stream in this 



area but our results agree at least qualitatively with their metal rich sample. 

The February 2008 sample of stars does not exhibit the same distinct signature we see 
in the August 2007 data. The velocities are much closer to disk velo cities and hav e a mu ch 



Yanny et al. 



(l2009f ) 



wider width. In Figure [13] we compare our RC results with those of 
There seems to be no division in the kinematics of the BHB stars based on abundance so all 
BHB stars are combined into a single histogram. While there is not a strong dynamically 
cold signal seen, there is a peak of stars in both samples with velocities > -100 km s~^. 
Both samples also have a tail extending toward more positive velocities. 



^ http : //model . obs-besancon . f r| 
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Shown in Figure [H] are histograms of GSR velocities in the area we probed with our 
August 2007 observat ions. Our RC results are again shown with the M giant r esults of 



Maiewski et al 



Law et al. 



(l2004f ) but the other curve is datsllj from the dynamical model of 
(j2005l ) which has been cut in location and distance to match our observations. One can see 
that there is not a large difference between the models. Our results tend to agree the best 
with the spherical model and agree the least with the prolate model (there is a large peak 
at positive velocities for stars with slightly larger distances in the prolate model and we see 
no trace of it in our data). Because the differences are small we cannot say with confidence 
that we definitely favor one model over another, however. It must be noted that the 



models of 



Law et al. 



( I2OO5I ) were constrained by the results of 



Majewski et al. 



(I2OO4I ) and 



therefore it is no surprise that they agree. It is significant, however, that our independent 



results n ot only match tho se of 



model of iLaw et al. 



Majewski et al. 



(12004) but a 



so ten d to rule out the prolate 



(I2OO9I )) suggests that our halo 



( 120051 ). A more recent paper (iLaw et al. 
may in fact be best characterized by a triaxial model instead of the axis ymmetric models 
used p reviously. A comparison with our data and the M giant results of 



Majewski et al. 



(I2OO4I ) is shown in Figure [T5i The same level of agreement can be seen in the triaxial and 



axisymmetric models in this area of the sky. 



the 



Figure flGl shows how our RC and the 



Law et al. 



Yanny et al. 



( l2009l ) BHB samples compare with 



( 120051 ) models constrained to match the area and distances of our February 
2008 observations. One can see here that while the two data sets agree, none of the models 
match the data. In particular, a strong positive peak is expected in all three halo models 
but is completel y absent from the data. When we compare this area of the sky to the 



triaxial model of 



Law et al. 



is similar to the prolate model of 



(I2OO9I) we see so r aethin g slightly different. The triaxial model 



Law et al. 



(l2005[ l with a couple of key differences. The 



^"available online at http: //www. astro. Virginia. edu/$\sim$srin4n/S gr7] 
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strong positive peak is still present but is shifted to lower velocities (slightly higher than 
km s~^). The negative peak is also shifted to smaller velocities. The result is a somewhat 
better agreement between our data, the BHB sample and the model. 

4.3. Density of RC Stream Stars 

As mentioned earlier, our sample is relatively complete with selection cuts only being 
made in magnitude and color. We therefore do not anticipate any kinematic bias in the 
RC sample. A small amount of metallicity bias is possible since the color of the RC is 
affected by metal abundance. However, the width of the color selection window is such that 
we expect this bias to be inconsequential. We therefore expect that we can make a fair 
determination of the RC density in the trailing and leading streams. 

Before determining this value we are also concerned with the question of whether the 
RC star sample is truly a sample of PopI, RC stars. This portion of the H-R diagram 
contains several stellar types which are not PopI stars. These include the metal poor 
K giants, the metal weak RHB stars and of course, potential contamination from the 
large population of main sequence stars. We have shown in section 14.11 that by removing 
the brighter stars in the sample the contamination of disk type stars is quite small. We 
further expect that the metal poor K giants are not a significant contaminate because the 
metallicity distribution has an average around [Fe/H] = -0.4. This is not indicative of K 
giants in the halo field. We further test this by plotting our August sample relative to the 
Besancon model, using the same color and magnitude cuts as our data and selected to 
match our observing fields. Figure [H] shows excellent agreement between our disk sample 
and that predicted by the Besancon model for a sample of stars chosen using our Tgjj/log g 
cuts. The peak velocity and dispersion are consistent, suggesting our disk sample is indeed 
composed of stars in the Galactic disk. The Besancon model halo stars, chosen from our 
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Te///log g cut, is quite different from the RC sample. The RC stars show a large peak 
with kinematics that are consistent with other publications of the trailing arm of the Sgr 
in this location. Furthermore, and most importantly, the RC sample shows virtually no 
population that would be consistent with the halo field. This indicates that the RC sample 
is constrained to the stream and is not represented in the general halo field. This would be 
expected if the RC sample is indeed predominately PopI, RC stars. Other than the primary 
stream peak, our sample only shows one very small peak, centered under the disk velocity 
distribution. We expect that this small peak may be a very small number of contaminating 
disk stars. 

With this information we are now ready to compute the density of RC stars in the 
trailing arm of the Sgr. We note that 3 of the 8 observed fields taken from the August run 
lie at the periphery of the stream as plotted in Figure [TJ These three fields also show no 
signal of the Sgr. In order to constrain the density as well as possible, we do not consider 
these three fields in the volume calculation and instead only concentrate on fields that are 
clearly located within the stream. The volume calculation is performed by summing up the 
total area observed in these five fields. Each WIYN field has a total area of 0.785 square 
degrees. This we multiplied by the 5 fields to get a total area of 3.95 square degrees. The 
depth of the observation was e stablished by using an e stimate of the appropriate absolute 



g magnitude from the Girardi (iGirardi fc Safaris 



20011) isochron es, and forcing the stream 



20091 ) for this region of the Sgr 



distance to match that of the new Law models (iLaw et al. 
stream. The optimal absolute g magnitude for the RC sample was found to be gabs = 0.6. 
Using this value we established that our stream data extended from 20 kpc to 38 kpc. 
The volume was computed by determining the volume of a cone truncated at 20 kpc. Our 
computed volume was V = 18.7 kpc^. 

The number of observed RC stars in the trailing arm of the stream is = 23 ± 3 
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stars. Because we were not able to observe every available star in the selected sample from 
the SDSS photometry, we used the ratio of the number of observed stars to that selected 
from the SDSS photometry as a filling factor. This gave us an expected number of Sgr stars 
of = 56 ± 6. This results in a density of p = 3.0 ± 0.3 RC stars/kpc^. 

Unlike the August data set, the February sample, which probes the leading arm of the 
Sgr, is surprisingly sparse. Figure [18] is the same as Figure [11] for the August data. It is 
readily apparent that a pronounced peak for the leading arm is absent in this figure. The 
disk stars of our sample again agree quite well with the location and dispersion predicted 
by the Besancon model data. The RC data, however appears to be scattered uniformly in 
velocity with only a very small peak near vgsr = -80 km s~^. This peak appears to be 



offset from the new 
is in contrast to the 



^aw model, 



Yanny et al 



'o r this location and distance, by about +30 km s ^. This 



( 120091 ) BHB data shown in Figure [13], which appears to 



sho w a very stro i ig sign al more closely matching the Law model. Both our RC sample and 



the 



Yanny et al. 



(|2009[ ) data follow virtually the same stripe across the Sgr stream. 



Compared to our August fields in the direction of the southern trailing arm, the 
northern leading arm has been an enigma. Although the total RC sample listed in Table 
[1] shows the February data to contain more RC stars than the August data, the February 
data shows a very weak signal and questionable peak kinematically. One possibility for the 
lack of a strong signal is proiection effects which increase the area of the stream on the sky. 



In the M giant sample of 



Majewski et al. 



fl2003[ ) the leading northern arm of the stream 
changes into what the authors call the northern 'fluff.' This is caused by a foreshortening 
effect that spreads the stream out as the stream descends on the local solar neighborhood. 
If foreshortening were to affect our data then the lack of a strong signal could be caused 
by the large area of sky over which the stream is spread as viewed from the Earth. This, 



however, does not seem to be the case. The region of the stream which we are sampling 
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is at a distance of 25 - 38 kpc with / ~ 24 5°. This is a region o 
associated with the northern arm from the 



Majewski et al 



the stream that is still 



( l2003l ) paper and a comparable 



distance to the Sgr sample in the trailing arm. 

To more quantitatively compare our August and February data, we have assumed that 
the small peak in February is actually a signal of the leading arm. The number of stars 
affiliated with the stream is = 12 ± 3. Using the selected SDSS photometry sample 
to assign a filling factor, the expected number of RC stars in the stream is = 24 ± 6. 
In order to sample the highest density of the stream we have used the six fields that are 
centered on the stream and computed a total area of 4.71 square degrees. The volume was 
computed in the same way as the August data, only over a distance of 25 - 38 kpc. This 
gave a total volume of 19.2 kpc^ and a density of p = 1.3 ± 0.3 RC stars/kpc^. Because 
it is not completely certain that the peak in February is associated with the leading arm, 
this density should be considered an upper limit. Even assuming the maximum value is 
correct, it is significantly less dense than the RC stars in the southern trailing arm. A table 
summarizing our star counts and densities can be found in Table [5l 

As a final test of the possible effects of foreshortening we have compared the density 
of the southern trailing arm and northern leading arm using the Law model. We chose 
areas from the model data that included our fields' coordinates and distance and computed 
the density in the same way as for the RC sample. The chosen fields were very uniformly 
populated by model points over the area in consideration. The ratio of densities between 
the southern trailing arm and the northern leading arm was found to be 1, suggesting that 
we should expect very similar densities in both arms. 

The scarcity of RC stars in the leading northern arm raises the questi on of whether 



there is a population gradient between the two streams. The fact that the 



Yanny et al 



(120091 ) data showed a strong signal in the PopII, BHB stars while we find a disparity 



- 20 - 



of RC stars seems to suggest that the northern leading arm is deficient in PopI stars. 



Majewski et al 



(|2003[ ) did not look for density variation s in the northern l eading arm 



because of the complications of foreshortening. However, 



Chou et al 



(120071 ) indentified a 



possible metallicity gradient in the M giants, with lower metallicities found for th e northern 



Bellazzini et al. 



fl2006h found 



arm when compared to the southern arm. Furthermore, 
evidence from photometric data that the BHB/RC ratio in the northern leading arm showed 
a clear decrease in RC stars compared to the Sgr proper. Their stream s ample was closer 



to the Sgr and in a location were the northern arm is well defined in the 



Majewski et al. 



( I2OO3I ) M giant sample (318 < / < 357). We expect that the effects of this gradient becomes 



even more severe as we move farther from the Sgr proper. 

The evidence for a gradient in the RC sample becom es even more dramatic when 



comparing to the new stream models by 



Law et al. 



fl2009h . Figures M and M show the 



predicted GSR velocity as a function of peri-galacticon passage (where 1 is the most recent 
passage) for the area of the stream surveyed in this paper. The solid vertical line is the 
mean velocity for our RC sample, while the dotted lines are the observed dispersion. Figure 
19] clearly shows that all the stars observed in the August data were deposited into the 
stream during the most recent passage of Sgr. This is in clear contrast to the February data 
shown in Figure [201 In this plot it is clear that the extra orbital distance from Sgr results 
in no stars being present from the most recent passage. There are, however, stream stars 
which have the correct GSR velocity that were deposited on earlier orbital passages. 

What does this say about a population gradient? It is important to real ize that the 
law m odels do not have any built-in population discriminate. The results of 



Yanny et al. 



20091 ) were for PopII, BHB stars, and as Figure [20] shows, it is well populated in the 
Northern leading arm. In contrast, our current RC sample has only a small suspected peak 
at ~ 80 km s~^. This small peak matches well with the velocity expected for stars deposited 
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two passages ago. Apart from this very weak signal, our data shows no other evidence for 
RC stars in the survey region of the leading arm. We do, however, find a very large RC 
signal in the August data. This strongly suggests that the young/intermediate, RC stars in 
the Sgr stream were preferentially stripped during the most recent peri-galacticon passage 
of the Sgr. 



5. Discussion 

We have shown that we can successfully separate RC type stars from a background 
of disk type stars using only stellar parameters obtained from medium resolution spectra 
and low to moderate signal to noise. This procedure applied to areas of the sky in the 
direction of the tidal streams of stars from the Sgr provides results that agree with previous 
measurements done on the stream. Most notably, the kinematics we have found match 
other populations of stars both qualitatively and quantitatively. Figure [2T] is a plot of 



Majewski et al. 



the G SR velocity as a function of the Sgr coordinate Asc as defined in 
( 2003 ) and transformed according to the resources available onlini^. Our August 2007 
sample is the group at Asc ~ 100° and the February 2008 sample is at Asc ~ 250°. The 
agreement between the three independent data sets is quite evident, as can also be seen in 



the histograms of Figures [12] and 

The metallicity distribution for our two samples also agree with previous measurements. 
We found values of [Fe/H] = -0.4 and [Fe/H] = -0.2 for the August 2007 and February 2008 
samples respectively. These match what has been found for the core of the Sgr and other 
dwarf galaxies. 

Since the RC population of stars in the Sgr system seems to be a sizeable sample. 



11 
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one should be able to probe differences in densities of stars stripped from the galaxy 
proper. Our comparison of the RC sample densities in the southern trailing arm and the 
northern leading arm indicates a factor of 2 disparity between the two stream locations. 
This evidence coupled with previous investigations suggests the young population stars 
have been stripped at a lower rate during previous orbits of the Sgr. This further suggests 
that the distribution of the young Sgr population was more centrally located than the 
older PoplI stars. Because our argument was based in part on stream models, it would 
be extremely beneficial to continue this type of program using complete samples of Popll 
stars, such as the BHB stars. This would allow a direct comparison of populations and 
firmly establish the significance of the population gradient. 

A clear benefit of studying the RC population in more detail is to further constrain the 
shape of the dark matter halo surrounding our Galaxy. When comparing to axisymmetric 
models, our data favors either a spherical or slightly oblate halo according to the 
August 2007 sample but we can not definitively distinguish between th e different shapes. 
Furthermore, our February 2008 RC sample and the BHB s ample from iYanny et al. 



Law et al. 



mm 



(120051 ). The agreement with a 



do not agree with any of the three models from 
triaxial model halo is similar to the spherical and slightly oblate axisymmetric models in 
the August 2007 sample. The February 2008 agreement is slightly improved when using 
a triaxial model but is still not nearly as well described as the August 2007 sample. We 
suspect, however, that the RC stars may be an important probe in helping to unravel the 
true shape of the halo. As seen in Figure EHl the BHB stars in the leading arm are likely 
composed of stars stripped over multiple passages of the Sgr, while the RC stars seem to 
only be found in the most recent passage, 1 and 2. This means that further observations 
of RC stars could help to constrain the most recent orbital passages of Sgr without the 
confounding contribution of multiple orbits present in the PopII sample. 
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The Sgr orbiting in and around the potential of our Galaxy offers a unique opportunity 
to study an isolated galactic system. The RC population can be used as an effective probe 
to help constrain both the shape of the dark matter halo and the initial distribution of stars 
in the Sgr galaxy. We have shown that this population of stars can be extracted from the 
significant background and yield results which are consistent with other populations of stars. 
Several different populations of stars with different ages are needed to fully characterize the 
interaction of the Sgr with our Galaxy and the RC is a justifiable and observable group to 
use as a young component. 

We would like to acknowledge the NOAO for granting time and providing funding to 
make the necessary observations for this research. 

Facilities: WIYN (Hydra). 
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Fig. 1. — Galactic coordinates for the spherical model of iLaw et al.l ( 120051 ) and the fields 
observed in this research. 
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Fig. 2. — Color-magnitude diagram of stars in our selected area of the sky. The box represents 
the selection criteria we used to select the RC program stars. 
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Fig. 3. — Histograms of the difference between three stellar parameters and the radial velocity 
of approximately 70 stars. The difference plotted is results obtained from the software 
developed for this research subtracted from high resolution results. 
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Fig. 4. — This shows the GSR velocity versus surface gravity for all program stars in the 
August 2007 observing area. 
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Fig. 5. — Surface gravity versus effective temperature plot for our program stars in the 
August 2007 observing area. The overlaid hues correspond to the cuts adopted to isolate 
stream stars from the sample. 
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Fig. 6. — Stars meeting the surface gravity and effective temperature requirements for the 
August 2007 observations are plotted here as GSR velocity versus the SDSS g magnitude. 
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Fig. 7. — The August 2007 stream subsample selected from the surface gravity and effective 
temperature cuts is plotted here as a function of the Galactic longitude. 
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Fig. 8. — Plot of GSR velocity versus surface gravity for the August 2007 stars which is 
similar to Figure HI Now the sample has been separated into disk (crosses) and stream 
(squares) stars using the surface gravity and effective temperature cuts. 
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Fig. 9. — Histogram of the metallicity for our August 2007 RC stars. The mean and char- 
acteristic width of the overlaid fit are given. 
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Fig. 10. — Histogram of the metallicity for our February 2008 RC stars. The mean and 
characteristic width of the overlaid fit are given. 
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Fig. 11. — Normalized histograms of GSR velocities for the August 2007 observations. The 
solid black histogram and curve are the data and fit respectively for stars selected as being 
Sgr stream RC stars. The black dashed histogram and curve are the data and fit respectively 
for stars determined to be disk type stars in the sample. The red dotted histogram and curve 
are the data and fit respectively for Besancon model disk stars in the same area of the sky as 
our observations. The blue dotted and dashed histogram and curve are for Besancon model 
halo stars. 
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Fig. 12. — Normalized histograms of GSR velocity for the are a probed by our 



observations. Our RC sa mple (black), t he M giant sample from 



and the BHB samples of lYanny et al. 



(!2nn"l the metal poor sample < -1.9 is blue and the 



Majewski et al 



Augus t 2007 



f l2004h (red), 



metal rich sample > -1.9 is green) are shown. 
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Fig. 13. — Normalized histograms of GSR velocity for the area probed by our Febr uary 2008 
observations. Our RC sample (black) and the BHB sample of lYanny et al.l (120091 ) (red) are 
shown. 
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Fig. 14. — Comparison of G SR velocities for our RC sample (solid black), t he M giant 



from 



Majewski et al. 



s ample 



Law et al. 



(120051 ) 



(120041 ) (dashed black), and the different halo models of 
(spherical is red, prolate is green and oblate is blue). These correspond to the area of our 
August 2007 observations. 
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Fig. 15. — Comparison of G SR velocities for our RC sample (solid black), the M gi ant sample 
from Majewski et al.l ( 120041 ) (dashed black), and the new triaxial halo model of iLaw et al. 
(120091 ) (red). This plot is for data from the August 2007 observations. 
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Fig. 16. — Compa r ison o f GSR velocities for our RC sample (solid black) , the BHB 



from 



Yanny et al. 



Law et al. 



s ample 



(120051 ) 



(120091 ) (dashed black), and the different halo models of 
(spherical is red, prolate is green and oblate is blue). These correspond to the area of our 
February 2008 observations. 
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Fig. 17. — Comparison of GSR velocities for our RC sample (solid black), the BHB sample 
from lYanny et al.l ( 120091 ) (dashed black), and the new triaxial halo model of lLaw et al.l (120091 ) 
(red). This plot is for data from the February 2008 observations. 
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Fig. 18. — Normalized histograms of GSR velocities for the February 2008 observations. The 
solid black histogram and curve are the data and fit respectively for stars selected as being 
Sgr stream RC stars. The black dashed histogram and curve are the data and fit respectively 
for stars determined to be disk type stars in the sample. The red dotted histogram and curve 
are the data and fit respectively for Besancon model disk stars in the same area of the sky as 
our observations. The blue dotted and dashed histogram and curve are for Besancon model 
halo stars. 
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Fig. 19. — GSR velocity as a function of peri-galacticon passa ge, where 
passage, for the August observing area. Points are from the iLaw et al. 
the sohd and dashed lines correspond to the mean and dispersion respectively for our RC 
sample. 
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Fig. 20. — GSR velocity as a function of peri-galacticon passa ge, where 1 i s the most recent 



passage, for the February observing area. Points are from the iLaw et ah 



mm model and 



the solid and dashed lines correspond to the mean and dispersion respectively for our RC 
sample. 
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Fig. 21. — GSR velocities are shown as a fu nction of the Sg r coor dinate Asc for our RC 
sample (solid boxe s), the M giant sana ple from lMajewski et al.l ( 120041 ) (solid circles), and the 
BHB sample from lYanny et al.l (120091 ) (open circles). 
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Table 1: Field Locations 





I 


b 


Field 


(degrees) 


(degrees) 


f02 


151.8 


-59.25 


f05 


148.9 


-59.64 


f07 


139.7 


-62.58 


f09 


157.3 


-57.44 


flO 


126.9 


-63.28 


fll 


164.0 


-56.68 


fl3 


105.1 


-61.06 


fl8 


171.7 


-51.74 


bOl 


249.3 


74.73 


b02 


246.3 


75.38 


b03 


252.0 


74.04 


b04 


243.1 


75.93 


b05 


254.5 


73.31 


b06 


239.7 


76.51 


blO 


231.7 


77.44 


bll 


260.9 


71.06 
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Table 2: Wavelength Regions and Weights 



Wavelength Range 








(A) 


Parameter 


Weight 


Notes 


4000 - 4080 


[Fe/H] 


1.0 




4500 - 4840 


[Fe/H] 


1.0 




4880 - 5000 


[Fe/H] 


3.0 




5200 - 5500 


[Fe/H] 


3.0 




4080 - 4120 


Te// 


1.0 


Hi 


4320 - 4360 


Te// 


1.0 


H-y 


4840 - 4880 


Te// 


2.0 




5000 - 5250 


logg 


1.0 


Mg Triplet 
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Ta ble 3: Comparison of Results 

High Resolution Difference SDSS Pipeline Difference 

Veiriable 





Mean 


Width (1 a) 


Mean 


Width (1 a) 


Te// (K) 


-15.71 


196.0 


17.53 


115.7 


Log g (dex) 


0.193 


0.265 


0.193 


0.356 


[Fe/H] (dex) 


-0.119 


0.240 


-0.059 


0.282 


Velocity (km s~^) 


5.721 


3.835 


7.859 


3.090 
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Table 4: Star Counts 



Data set Aug 07 Feb 08 

Disk 263 225 

Stream 43 58 

Total Observed 306 283 

SDSS Total 751 602 
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Table 5: RC Star Counts and Derived Densities 







RC Density 


Area 


RC Number 


(stars/kpc^) 


Aug 07 


56 ± 6 


3.0 ± 0.3 


Feb 08 


24 ± 6 


1.3 ± 0.3 



